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ABSTRACT 

We describe a near-infrared imaging survey covering ~ 1 deg^ of the 
Chamaeleon I dark cloud. The survey is complete for K < 15.0, H < 16.0, 
and J < 16.5, roughly two magnitudes more sensitive than previous large scale 
surveys. We use the large number of background stars detected to derive an accu- 
rate near-infrared extinction law for the cloud and select new candidate members 
with near-infrared color excesses. We list ~ 100 candidates of the cloud with K 
> 12.0, based on their positions in the J — H, H — K color-color diagram. These 
new stars have low luminosities (K 12 - 16, H — K > 0.5 - 1.5) and may have 
masses close to or even below the hydrogen burning limit. 

Subject headings: ISM: individual (Chamaeleon I) — ISM: dust, extinction 
— star: formation — stars: low-mass stars, brown dwarfs — star: pre-main 
sequence 

1. Introduction 

The Chamaeleon I dark cloud (a ~ 11*^, 5 ~ —77°) is an active stellar nursery with 150 
or more known young stars (Schwartz 1991; Gauvin and Strom 1992; Lawson et al. 1996; 
Cambresy et al. 1998; Comeron et al. 1999). Located at a distance of ~ 160 pc (Whittet et 
al. 1997), it has a relatively small angular size (Boulanger et al. 1998, ~ 3 deg^). Because 
of its proximity and moderate extension on the sky the Chamaeleon 1 star-forming region is 
an attractive place to attempt to measure the initial mass function (IMF), particularly at 
low masses (i.e., for masses < 0.1 Mq). 

Here we report results of a JHK imaging survey of the Cha I dark cloud, which covers 
an area of ~ 1 deg^ and is complete for K < 15.0, H < 16.0, and J < 16.5. Our survey 
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complements the IJKs observations of DENIS (Deep Near Infrared Southern Sky Survey, 
ESO, Im tel, Epchtein (1997)) described by Cambresy et al. (1998). The DENIS survey 
covers a larger area, 1.5° x 3°, of the cloud; but, with a sensitivity of Kg < 13.5, it is 
less sensitive to the lowest mass cloud members. Cambresy et al. (1998) proposed 54 new 
candidate young stars, based on their locations in the I — J / J — K diagram. In addition our 
survey area comprises two areas, in the the northern region, recently observed by Persi et 
al. (1999) and Oasa et al. (1999). The data reported by these authors are ~ 2 magnitudes 
deeper than ours but cover modest areas on the cloud (~ 3' x 3', Persi et al. (1999); ~ 6' x 
6', Oasa et al. (1999)) centered close to the high-velocity C^^O bipolar outflow previously 
detected by Mattila (1989). Persi et al. (1999) identified a new Class I source of the cloud. 
Oasa et al. (1999) proposed 9 new very low luminosity members of the cloud. 

Our survey recovers 40 of these DENIS candidates (lying on common areas to both 
surveys) and the stars with K < 15.5 detected by Persi et al. (1999) and Oasa et al. 
(1999). In addition the present observation yields ~ 100 new candidates with K > 12 and 
near-infrared excess emission in the J — H,H — K color-color diagram. The low luminosities 
of these candidates suggests they may be very low mass cloud members, with masses close 
to or below the hydrogen burning limit. 

We describe our observations, data reduction, and analysis in §2. The astrometry and 
photometry were obtained for ~ 11,090 sources at JHK detected in our survey region. In 
§3 we derive a reliable extinction curve for background stars and use this curve to propose 
~ 100 new potential young stellar objects based on their location in the J — H / H — K 
diagram in §4. We conclude with a brief summary in §5. 
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2. Observations and Data Analysis 

We obtained JHK imaging data for the Chamaeleon I dark cloud and three relatively 
unreddened control fields on 13-16 February 1995, 8-11 March 1996, 16-19 April 1997, and 
1-2 November 1998 with CIRIM (the Cerro Tololo Infrared Imager) at the CTIO 1.5m 
telescope. The CIRIM uses a 256 x 256 HgCdTe NICMOS 3 array which provides a field 
of ~ 4.9' X 4.9' with a plate scale of 1.16" per pixel. We covered an area of ~ 0.65° x 1.5° 
on the cloud on a rectangular grid, with 1' overlap between adjacent frames. We acquired 
two 6x5 sec exposures for each field, shifted by 20". Figure 1 shows the extent of the 
survey relative to the H2CO contour map measured by Toriseva and Mattila (1985) and the 
distribution of previously known young members of the cloud (Schwartz 1991; Gauvin and 
Strom 1992; Lawson et al. 1996; Cambresy et al. 1998; Comeron et al. 1999). Cha I has a 
remarkable low-mass star formation activity localized mainly in three groups. These small 
clusters are situated near three refiection nebulae: Ced 112 (HD 97300), Ced 110, and Ced 
111 (HD 97048), from North to South. 

Our three control fields, covering an area of ~ 0.05 deg^ and selected from a visual 
inspection of the ESO Red Sky Survey prints, are close to the survey region and relatively 
free from significant optical extinction. These off-cloud regions lie outside the limits of 
Figure 1, ~ 1.5-2 deg away from the H2CO boundaries, one to the North and the other two 
to the East and West, respectively. Table 1 gives positions and total areas covered by each 
of the off-cloud regions. 

We processed the data using standard techniques with the software package IRAF^. 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc. under contract to the 
National Science Foundation. 
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Prom each program frame, we subtracted an average dark frame, divided by a normalized 
fiat-field frame appropriate for each filter, and then subtracted a flattened sky frame. The 
dark frame is the average of ~ 40 individual dark images taken at the beginning and end 
of each night. We constructed flat- field frames for each night by median-filtering all 6 x 5 
sec frames taken in each filter. We also obtained a set of " dome flats " in each filter before 
observing each night using an illuminated white spot in the CTIO 1.5-m dome. Flat-fields 
constructed from " dome flats " and " sky flats " were identical. Sky frames were generated 
from a median-filtered set of 20-30 fiattened frames in each filter. In this case, we combined 
individual object frames obtained close in time and in position to the program frame. 

To produce a combined image from each dithered pair of images, we aligned the frames 
using the IRAF subroutines GEOMAP and GEOTRAN, added the co-aligned frames 
together, and trimmed the resulting image to remove bad pixels at the edges of the frames. 
We selected the K images as the reference frames and transformed the J and H images 
to the same pixel scale as the K images. We used DAOFIND to locate stars Aa above the 
local background and added to the DAOFIND list all stellar objects missed by this routine, 
found by visual inspection of each image. We then derived photometry for each image using 
the APPHOT PHOT task, using a circular aperture with 5" radius. This aperture size 
includes the total fiux for the large majority of sources. However, in some cases, we used a 
smaller aperture (i.e. 2.5") to avoid contamination from very close stars. Several iterations 
of this process produced an homogeneous set of data for the complete survey. We detected 
~ 11,090 sources at JHK in our survey region. 

To cahbrate our photometry, we observed on each night a set of 10-15 standards from 
Elias et al. (1982) and from the UKIRT faint JHK standard stars hst (Casali and Hawarden 
1992). We estimate an uncertainty of ± 0.03 mag in our calibration. The standard stars 
were observed at a similar airmass range as our target fields. Airmass corrections for our 
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data are smaller than our photometric uncertainties. We used common stars on dithered 
pairs and measurements of duplicate stars in the overlapping regions of adjacent frames to 
estimate photometric uncertainties for the stars in our survey region. Table 2 lists typical 
photometric uncertainties in each magnitude bin. These differences reach 0.3 mag for J ~ 
17, H ~ 16, and K ~ 15. The 5(7 limiting magnitudes are K — 14.5, H — 15.5, and J = 



To make a check on the completeness of the survey, we constructed a hst of sources 
from the three off-cloud regions (see Table 1). Figure 2 shows the magnitude histograms of 
background stars as function of the JHK magnitudes. The number of stars per bin increases 
monotonically up to the completeness limit and then turns over. The corresponding limits 
in each filter are: 16.5 at J, 16.0 at H, and 15.0 at K. 

To compare our photometry with results for previous surveys, we observed ~ 50 stars 
projected on Bok globule 2 in the Coalsack (Tapia 1975). Jones et al. (1980) and Elias et 
al. (1983) previously reported photometry of these stars and derived a comparison between 
the Mt. Stromlo/AAO photometry of Jones et al. (1980) and the CIT photometric system 
of Elias et al. (1982). The common photometric data between Jones et al. (1980) and our 
survey, together with the transformations listed in Elias et al. (1983), allow us to derive a 
transformation between the CIRIM and CIT photometric systems. Using the Press et al. 
(1992) routine FITEXY, straight hue fits to the common data yield: 



16.0. 



KciRiM = +0.03 ± 0.02 + 0.99 ± 0.01 • Kcit 



(1) 



(J — H) CIRIM 



0.04 ±0.03 + 0.95±0.02-(J-H)ciT 



(2) 
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(H-K)ciRiM = -0.02 ±0.02 + 1.01 ± 0.02 • (H-K)ciT (3) 

The quoted uncertainties in tliese relations are the la errors from our fits to the Mt. 
Stromlo/AAO data and those from the Ehas et al. (1983) fits, added in quadrature. The 
differences between the CIRIM and CIT system are neghgible for K and H — K, in 
agreement with results quoted in the CIRIM manual written by R. Elston and J. Ehas^. 
We thus assume that the natural CIRIM system is identical to the CIT system for K and 
H — K. Our derived color term ior J — H agrees with Elston & Ehas; the CIRIM J — H 
color is bluer than the CIT J — H due to its different lens design. We apply the color term 
to our J — H data and thus quote colors in the CIT system. 

To derive coordinates for our survey stars, we used WCSTool^ (Mink 1997), a suite of 
programs to calculate a direct transformation between the coordinates of the image (x,y) 
and the sky coordinates (a, 5). We measured transformation coefficients adopting matches 
between program stars and stars in the U.S. Naval Observatory SAl.O Catalogue^. These 
matches provided a good transformation for ~ 90% of the frames. For the rest of the 
frames, usually corresponding to the most obscured regions, we found few or no matches 
between our frames and stars in the catalog. In these cases, we used stars in common with 
adjacent frames with good coordinates to obtain coordinates for stars on unmatched frames. 
When an individual source was detected in more than one filter we adopted the average 
coordinate. 

^Available at http:/ /www. ctio.noao.edu/instruments/ir_instruments/cirim/cirim.html 
^Available at ftp://cfa-ftp.harvard.edu/pub/gsc/WCSTools 

^The U.S. Naval Observatory SAl.O Catalogue is distributed by the U.S. Naval Observa- 
tory, Washington DC. 
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To verify the coordinates, we measured coordinate differences between stars on 
adjacent frames. Many stars have two good coordinates; others have as many as four. After 
visually inspecting the duplicate stars and rejecting stars too close to the edges of frames, 
we averaged coordinates of the remaining duplicates and derived la errors. Finally, we 
compared coordinates derived from the Digitized Sky Survey^ (DSS) plates. We estimate 
an average uncertainty in our positions of 1". These errors are as large as 2" for sources in 
the most obscured regions of the cloud. 

Having verified the quality of the coordinates and photometry for our survey, we 
search for pre-main sequence stars with near-infrared excess emission in the cloud using the 
near-infrared color-color diagram. To identify these stars, we derive a reliable extinction 
curve for background stars in §3 and then consider the near-infrared color-color diagram in 
§4. 

^Based on photographic data obtained using The UK Schmidt Telescope. The UK 
Schmidt Telescope was operated by the Royal Observatory Edinburgh, with funding from 
the UK Science and Engineering Research Council, until 1988 June, and thereafter by the 
Anglo-Australian Observatory. Original plate material is copyright (c) the Royal Observa- 
tory Edinburgh and the Anglo- Australian Observatory. The plates were processed into the 
present compressed digital form with their permission. The Digitized Sky Survey was pro- 
duced at the Space Telescope Science Institute under US Government grant NAG W-2166. 
Copyright (c) 1993, 1994, Association of Universities of Research in Astronomy, Inc. All 
right reserved. 
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3. The Cha I Near-Infrcired Reddening Law 

To derive the near-infrared reddening law for the Cha I dark cloud, we follow Kenyon 
et al. (1998a), who developed a generalized photometric technique to compare the colors 
of reddened stars with the colors of nearby, " unreddened " comparison stars. Kenyon et 
al. (1998a) assume that the stellar population behind the cloud is identical to the stellar 
population in off-fields several degrees away. They derive the J — H and H — K color 
excesses for each on-cloud source relative to every off-cloud source and then compute the 
average and median color excesses for each on-cloud source. The probable error of the 
average color excess is the sum in quadrature of the errors of the on-cloud and off-cloud 
colors. For the median color excess, the probable error is the inter-quartile range. We 
divide the Cha I sources into a Complete Sample containing all stars with K < lA and a 
Restricted Sample, where we remove known young stars in the cloud from the Complete 
Sample. The sample of known young stars includes 94 sources (lying on our survey region) 
from the literature (Schwartz 1991; Gauvin and Strom 1992; Lawson et al. 1996; Cambresy 
et al. 1998; Comeron et al. 1999). This sample probably is not complete; however, any 
incompleteness does not affect our analysis significantly. 

Figure 3 shows average color excesses of Cha I stars for two K < 14 samples, the 
Complete Sample and the Restricted Sample. The colors have been transformed to the 
CIT system using equations 1-3. The color excesses are highly correlated, with a Spearman 
rank-order correlation coefficient of = 0.65 for both samples. The probability for no 
correlation between the two color excesses is formally zero according to the Spearman 
rank-order test. Straight line fits to color excess measurements with Eh-k < 2.0 yield 
Ej_h/Eh-k = 1-63 ± 0.03 for the Complete Sample and Ej_h/Eh-k = 1-76 ± 0.02 for 
the Restricted Sample. We derive identical slopes using median color excesses for each 
source instead of the average color excesses. 
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Our extinction results for the Restricted Sample are identical to the Bessell and 
Brett (1988) near- infrared extinction law, Ej_h/Eh^k = 1-75 transformed to the CIT 
photometric system. The Cha I results differ from the near-infrared extinction law for p 
Oph, Ej^h/Eh-k = 1-57 ± 0.03 (Kenyon et al. 1998a; Ehas 1978, for example) at more 
than the 3a level, and are inconsistent with results from the He et al. (1995) survey of 
luminous southern stars, Ej_h/ E^-k = 1-47 ± 0.06, at roughly the Scr level. 

Results derived for the Complete Sample demonstrate that pre-main sequence stars 
with near-infrared excesses skew the extinction law to smaller values (see also Kenyon et 
al. (1998a) and references therein). Figure 3 shows many previously identified sources with 
apparent near-infrared excesses. These sources probably skew our results for the reddening 
law derived from the Restricted Sample, although this effect should be small. Kenyon et al. 
(1998a) described a method, using the " reddening probability distribution ", to correct the 
reddening law for previously unidentified pre-main sequence stars and to derive a sample of 
candidate pre-main sequence stars with near-infrared excess emission. We now apply this 
method to the Restricted Sample of Cha I sources. 

Kenyon et al. (1998a) defined the reddening probability distribution, p(E^,Ey), as the 
probability of measuring a pair of color excesses, Er^ and Ey, where x and y are color indices 
and N — J p{Ex, Ey) dE^dEy is the number of reddening measurements. The number, N, 
is also the number of on- field stars in the sample. The density function, p{Ex, Ey), depends 
on the distributions of colors in the on- field and the off-field. If poff{x,y) is the off-field 
color distribution and Pon{x,y) as the on-field color distribution, the reddening probability 
distribution is: 



where — i + ki — io and k^ — j + k2 — jo- The integers io and Jq measure the zero point 




(4) 



fci=l k2=k, 



offsets of tfie color grids; indices i and j span the full range of color indices x and y. The 
color distributions use a the kernel density estimator: 



(5) 



where h is the smoothing parameter and K is the kernel (Silverman 1986). We follow 
Kenyon et al. (1998a) and adopt 



for fast computation (x = (,x, Xi, y, yi)). We adopt h = 0.2, which is roughly twice the 
la error of our photometry at the K — 14 mag survey limit. Smaller values for h produce 
noisy grids; larger values are inconsistent with the photometric errors. 

Figure 4 shows the color distribution functions for the off-field and the Complete 
Sample (CS). The off- field density in the left panel has a weak maximum at the color 
expected for G-type giants and a sharp maximum at the color expected for K-type dwarfs. 
There is a plateau consistent with colors for M-type dwarfs and an extension of the dwarf 
locus towards the near-infrared colors of K- and M-type giant stars. The Cha I color 
distribution in the right panel peaks close to the off-field maximum but is much more 
extended. The long axis of the contours roughly follows the reddening line derived above, 
and there is a weak extension of sources with near-infrared excess emission at H — K — 
0.75 and J - = 1.75. 

Figure 5 shows the reddening probability functions constructed from the color 
distributions for the Complete Sample (CS) and the Restricted Sample (RS). We derive the 
slope of the reddening law which follows the major axis of each contour in these diagrams 




otherwise 



(6) 
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using the technique outhned in (Kenyon et al. 1998a, equations 12-15). This analysis yields 
Ej_h/Eh^k = 1.77 ± 0.03 for the Complete Sample and Ej_h/Eh^k = 1.80 ± 0.03 
for the Restricted Sample. These results are consistent with the reddening law derived in 
Figure 3 given the la errors of each fit. 

4. New Candidate Pre-Main Sequence Stars 

Figure 6 shows the H — K / J — H diagram for our complete survey region. As in 
Figures 3-5, we display stars with K < 14.0, where the photometric errors are < 0.1 mag. 
The solid line corresponds to the main sequence locus (Bessell and Brett 1988); the dashed 
lines, parallel to the reddening vector derived in §3, define the reddening band extending 
from the main sequence. 

We distinguish two groups of sources in Figure 6: a) stars that he along or follow the 
reddening band direction, and b) sources that show near-infrared excess, located to the 
right of the reddening band. The first group, comprising the majority of the sources, show 
negligible near-infrared color excess, and can be dereddened to lie close to the standard 
dwarf sequence. Most of these sources are probably behind the cloud^. The color excess 
of the second sample cannot be attributed to the reddening by dusty material in the 
cloud. Several authors have found that circumstellar material (disks -|- inf ailing envelopes) 
surrounding the young central star can produce these color excess at the near-infrared 
wavelengths (Lada and Adams 1992; Kenyon and Hartmann 1995). Formally speaking 
we have detected ~ 300 sources lying to the right of the reddening band in our survey 
with K < 14.0. Roughly 1/3 are known members of the cloud and DENIS candidates from 

^The WTTS that belong to the cloud and that, as a class of young stellar objects, has 
no significant near-infrared excess may be included in this group. 



-13- 



Cambresy et al. (1998). The other 2/3 are new detections. In Figure 6 we have used starred 
symbols to indicate new candidates with significant near- infrared excess (i.e., H — K > 
0.5). This group comprises ~ 50 objects. Table 3 gives coordinates and magnitudes for 
these stars. The rest of the sources either lie very close to the reddening line or to the bulge 
defined by the background stars. 

Finally, our survey has also produced a list of ~ 50 faint objects (i.e., K > 14 - 16 ), 
with large near- infrared excess (i.e., H — K > 0.8). Figure 7 shows the location of these 
stars in the color-color diagram. They are located well to the right of the main sequence 
locus. In this manner even allowing typical photometric uncertainties for these magnitudes 
(i.e., ~ 0.3 mag for K ~ 15) the sources still have color excesses. Table 4 gives coordinates 
and magnitudes for these stars. For these faint sources, we checked each of the images 
visually to avoid possible confusion with cosmic ray events, bad pixels, plate artifacts, or 
contamination with close stars. 

A few of the objects listed in the Table 3 and Table 4 were only detected in two of the 
filters {K and H) and therefore J must be > 18.0. These are extremely red objects and 
thus good candidate members of the cloud. Some of the objects display extended or fuzzy 
images principally at K. Four candidates were also detected by the ISOCAM survey of 
the Chamaeleon I dark cloud and proposed as new low-mass candidate members based on 
their position on the mid- infrared color-color diagram (Persi et al. 2000). We also have in 
common two stars with Oasa et al. (1999) (see Table 3 and Table 4). 

Sources with near-infrared excess emission are either young stars in the cloud, galactic 
sources such as planetary nebulae (e.g., Gughelmo et al. (1998)), or galaxies. Extragalactic 
source counts predict 30-60 galaxies in an area similar in size to our survey region for K < 
14 (e.g., Szokoly et al. (1998); Vaisanen et al. (2000) and the references therein). The 
near-infrared extinction towards Cha I is significant (Cambresy et al. 1997); we thus expect 
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~ 30% to 50% fewer galaxies in our survey compared to unobscured regions. Galaxies also 
tend to have smaller near-infrared excesses compared to young stars. Most galaxies have 
colors similar to red main sequence stars, J — H ^ 0.6-0.8 and H — K 0.25-0.45, with 
modest near-infrared excesses of 0.2-0.3 mag in H-K (e.g., Progel (1985); Impey et al. 
(1986)). Unpublished 2MASS data support this conclusion for nearby galaxies (see Jarrett 
et al. (1998)^). Our survey is not deep enough to detect more distant galaxies which could 
have redder near-infrared colors than measured for the nearby galaxy samples. 

Several types of galactic sources, such as carbon stars and planetary nebulae, can 
also lie outside the reddening band in the near-infrared color-color diagram. The surface 
densities of these objects towards Cha I should be small, < 5-10 deg~^ for carbon stars and 
< 1-2 deg~^ for planetary nebulae (e.g., Gughelmo et al. (1998); Ortiz and Maciel (1994)). 

The spatial distribution of the near-infrared excess sources provides the best evidence 
that these objects are pre-main sequence stars. Figure 8 compares the spatial distribution 
of the new candidate objects with the positions of previously known young stars (Schwartz 
1991; Gauvin and Strom 1992; Lawson et al. 1996; Cambresy et al. 1998; Comeron et al. 
1999) and the H2CO contour maps (Toriseva and Mattila 1985). The new candidates 
tend to cluster around Ced 112 (HD 97300), Ccd 110, and Ced 111 (HD 97048) in a 
manner similar to the previously known members of the cloud (see Figure 1). We expect 
extragalactic and other galactic sources such as carbon stars and planetary nebulae to avoid 
regions of high obscuration and see no evidence for this tendency in our candidates. 

Assuming an average inter-cloud extinction of Ay ~ 5 (Cambresy et al. 1997; Cambresy 
1999, i.e., Ak ~ 0.5), a main-sequence star at the H burning limit and at the distance 
of the Cha I cloud (Whittet et al. 1997, 160 pc), would have K - 15.5 - 16.0 (Henry 

^Available also at http://spider.ipac.caltech.edu/staff/tchester/2mass/analysis/galaxies/colors/ 
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and McCarthy 1993). Low mass members of the Chamaeleon I dark cloud must be much 
younger, still contracting toward the main sequence, and thus significantly more luminous 
than main sequence stars of 0.1 Mq. Our imaging survey is complete to K ~ 15.0 and 
sensitive down to 16.5. The objects detected in this survey (with K 12 - 16, H — K > 
0.5 - 1.5) are among the lowest luminosity and, presumably, lowest mass members of the 
cloud. Some of them are near or, probably, below the H burning limit. 

Some of our candidates (with X ~ 12 - 14) have optical counterparts on the DSS 
plates. We crudely estimate V ~ 17-22 for this subgroup of candidates. Moderate 
resolution optical spectra would allow a search for Li I 6707 absorption, an indicator of 
youth, as well as other atomic (such as Ha) and many forbidden emission lines usually 
present in the spectra of young stellar objects (Kenyon et al. 1998b, for example). To 
investigate further the nature of the rest of our candidate objects we need additional 10 /xm 
photometry and/or near- infrared spectroscopic observations. In particular, near- infrared 
spectra would allow us to place the stars in the HR diagram and thus, adopting recently 
developed pre-main sequence tracks (D'Antona and Mazzitelli 1997; Palla and Sthaler 1999, 
for example), determine masses and ages for the individual sources. Mass determinations 
are needed to estimate the lower end of the IMF. Present determinations of the IMF in 
the Cha I cloud are hmited to 0.3 Mq and for masses < 0.6 Mq these determinations 
are uncertain (Appenzeller et al. 1983; Lawson et al. 1996). Age estimations will provide 
additional support to the concept of a coeval star-formation process in this cloud (Gauvin 
and Strom 1992). 

5. Summciry 

We have used the large number of background stars detected in our survey region to 
derived a reliable near-infrared extinction law for the Chamaeleon I cloud. Our analysis 
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yields Ej^h/Eh-k — 1-80 ± 0.03 for the Restricted Sample (i.e., eliminating previously 
known pre-main sequences stars). This result differs with the Kenyon et al. (1998a) 's 
determination for the p Oph cloud {Ej_h/ Eh-k — 1-57 ± 0.03) at more than 3cr level. 
Based on preliminary results for other two cloud (Taurus and IC 348), in addition to 
their determination for p Oph, Kenyon et al. (1998a) suggest that a real variation in the 
near-infrared reddening law occurs from region to region. Our results for Cha I provide 
additional support for this variation. 

The ratio of total to selective extinction, Ry, seems also to change for different clouds 
(see Kenyon et al. (1998a) and the references therein). In addition, Ry varies across the 
Cha I dark cloud (from Ry ~ 3 to Ry ~ 5.5; Whittet et al. (1997), see also Hyland et al. 
(1982)). For the p Oph, Vrba et al. (1993) obtained a mean value of Ry ~ 4, although 
individual determinations vary within a range similar to that found for the Cha I dark 
cloud. Hayakawa et al. (1999) recently found that the Ay - N{CO) relation also varies 
among different star- forming clouds and suggest that Av/N{CO) may also change across 
the Chamaeleon I dark cloud. These results indicate different dust-grain properties, in 
particular grain-sizes, for different nearby star-forming regions and caution against using the 
" universal " reddening law to derive extinction estimates for young stars in star-forming 
regions. 

Previous near-infrared observations of the cloud focussed on IRAS-selected sources 
(Assendorp et al. 1990; Whittet et al. 1991; Prusti et al. 1991) or already known and 
candidate young stellar objects detected at optical and X-ray wavelengths (Gauvin and 
Strom 1992; Hartigan 1993; Lawson et al. 1996; Comeron et al. 1999). These investigations 
were sensitive to objects with K ~ 11 and limited to small areas around the target objects. 
Present imaging near-infrared surveys of the cloud have overcome these limitations (both 
in sensitivity and spatial coverage) and greatly increased the number of potential members 
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of the cloud. (Persi et al. 1999) detected one new Class I source of the cloud with K ~ 13. 
Oasa et al. (1999) proposed 9 new candidates with K = 13 - 16. (Cambresy et al. 1998) and 
this paper carry out large scale near-infrared survey on the cloud. (Cambresy et al. 1998) 
proposed 54 new candidates members with K < 13. This paper presents an additional set 
of ~ 100 new stars with K ~ 12 - 16. In addition to these candidates selected on basis of 
their near-infrared color excesses, Persi et al. (2000) have recently proposed a list of 74 new 
low mass candidates detected by ISOCAM and selected on basis of their color excesses at 
the mid-infrared wavelengths. 

These objects are among the lowest luminosity and, presumably, lowest mass members 
of the cloud. Some of them are near or, probably, below the H burning limit. A spectroscopic 
follow up of the optically visible candidates should provide additional indications of their 
pre-main sequence nature. Optical and/or near-infrared spectral types of our proposed 
candidates, in combination with recently developed pre-main sequence evolutionary tracks, 
will allow us to determine masses and ages. Reliable mass determinations are required 
to calculate the IMF, particular towards the sub-solar and probably into the sub-stellar 
mass regime. Ages determination will help to reconstruct the history of the star-formation 
process of the cloud. 
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and the National Aeronautics and Space Administration (grants: NAGW-2919 and 
GO-06132.01.94A). M.G. acknowledges support from the National Science Foundation 
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Fig. 1. — The extent of our near-infrared survey relative to the H2CO contour maps from 
(Toriseva and Mattila 1985) and the distribution of previously known members of the cloud 
(Schwartz 1991; Gauvin and Strom 1992; Lawson et al. 1996; Cambresy et al. 1998; Comeron 
et al. 1999, dots). The large starred symbols indicate the positions of Ced 112 (HD 97300), 
Ced 110, and Ced 111 (HD 97048) from North to South, respectively. 

Fig. 2. — Histogram distribution of magnitudes for background stars in our control fields. 
The completeness limit in each band is indicated by the dashed vertical lines. 

Fig. 3. — Average color-excess diagram for two samples of the Chamaeleon I stars with K 
< 14.0. The left panel shows the Complete Sample (CS) and the right panel the Restricted 
Sample (RS). This second sample only includes near- infrared sources not known to be pre- 
main sequence stars. 

Fig. 4. — Color density distribution for sources in the off-field region (left panel) and for 
Cha I sources Complete Sample (CS, right panel). 

Fig. 5. — Reddening probability distribution for the Complete Sample (CS) and the Re- 
stricted Sample (RS). The dot-dashed hne indicates our best reddening law; the dashed hne 
shows a standard reddening law with Ej_h/Eh-k = 1-75. 

Fig. 6. — Color-color diagram for near-infrared sources with K < 14.0 in our survey region 
(see Figure 1). The solid line indicates the main sequence locus (Bessell and Brett 1988) 
and the two parallel lines define the reddening band, extending from the main sequence. To 
define this band we used the reddening vector {Ej^h/Eh-k = 1-80 ± 0.03) derived for the 
Restricted Sample (RS in Figure 5, left panel). Typical photometric errors for K ~ 14.0 are 
displayed in the upper left corner. 

Fig. 7. — Color-color diagram for near- infrared sources with K > 14.0 in our survey region 
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(see Figure 1). The solid line indicates the main sequence locus (Bessell and Brett 1988) 
and the two parallel lines define the reddening band, extending from the main sequence. To 
define this band we used the reddening vector {Ej_h/Eh-k = 1.80 ± 0.03) derived for the 
Restricted Sample (RS in Figure 5, left panel). Typical photometric errors for K ~ 15.0 are 
displayed in the upper left corner. 

Fig. 8. — Spatial Distribution of the new candidate young stellar objects (triangles) fisted in 
the Table 3 and Table 4 in relation to the H2CO contour maps (Toriseva and Mattila 1985) 
and the previously known members of the cloud (Schwartz 1991; Gauvin and Strom 1992; 
Lawson et al. 1996; Cambresy et al. 1998; Comeron et al. 1999, dots). The large starred 
symbols indicate the position of Ced 112 (HD 97300), Ced 110, and Ced 111 (HD 97048) 
from North to South, respectively. 
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Table 1. Off-Cloud Regions. 



Region a(2000.0) 5(2000.0) Area 



I 11 31 43 -77 38 22 5'x5' 

II 10 41 24 -77 38 35 lO'xlO' 

III 11 16 52 -77 08 42 5'xlO' 



Note. — Units of right ascension are hours, 
minutes, and seconds, and units of declination 
are degrees, arcminutes, and arcseconds. 
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Table 2. Mean Photometric Errors. 



Mag. Bin K H J 



8 - 9 0.02 0.01 0.01 

9 - 10 0.03 0.02 0.02 

10 - 11 0.03 0.02 0.02 

11 - 12 0.03 0.03 0.03 

12 - 13 0.04 0.03 0.03 

13 - 14 0.07 0.04 0.03 

14 - 15 0.12 0.08 0.05 

15 - 16 0.26 0.12 0.07 

16 - 17 0.46 0.20 0.09 
17-18 ••• 0.38 0.15 
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Tablc 3. Candidate Pro-Main Scqiicncc Stars with K < 14.0. 



Star q;(2000.0) 5(2000.0) K H-K J-H Id. Ref. 



1 


11 


04 34.0 


-78 


01 


26 


13.34 


0.68 


1.00 


2 


11 


04 41.5 


-77 


38 


52 


13.64 


0.73 


1.09 


3 


11 


04 49.2 


-78 


04 


46 


12.48 


0.75 


1.03 


4 


11 


04 59.5 


-77 


17 


00 


12.59 


0.85 


0.84 


5 


11 


05 27.2 


-77 


33 


09 


11.80 


2.99 


1.24 


6 


11 


05 43.0 


-77 


29 


13 


13.79 


1.34 


1.74 


7 


11 


05 43.4 


-77 


21 


22 


13.81 


0.83 


0.85 


8 


11 


05 48.8 


-76 


40 


17 


13.94 


0.73 


0.55 


9 


11 


05 57.8 


-76 


40 


56 


13.06 


0.84 


0.79 


10 


11 


06 08.8 


-76 


56 


02 


13.56 


0.73 


0.91 


11 


11 


06 18.6 


-77 


57 


53 


13.95 


0.89 


1.42 


12 


11 


06 27.7 


-77 


35 


54 


13.26 


1.10 


1.65 


13 


11 


06 39.3 


-77 


36 


04 


12.40 


1.23 


1.78 


14 


11 


06 40.1 


-76 


48 


35 


12.48 


0.94 


1.03 


15 


11 


07 06.5 


-76 


37 


17 


10.22 


1.50 


0.93 


16 


11 


07 18.8 


-77 


31 


56 


12.63 


1.93 


1.49 


17 


11 


07 27.9 


-76 


48 


35 


13.91 


0.74 


0.38 


18 


11 


07 46.9 


-76 


15 


17 


12.22 


0.71 


0.94 


19 


11 


07 50.0 


-77 


38 


07 


13.50 


1.55 


2.23 


20 


11 


08 20.6 


-76 


26 


49 


13.80 


0.67 


0.62 



ISO-Chal 79 1 
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Table 3 — Continued 



Star 


o(2000.0) 


^K2000.0) 


K 


H - K 


J - H 


Id. 


Rcf. 


21 


11 


08 22. D 


-76 49 19 


12.32 


1.12 


0.82 






22 


1 1 
11 


no /II c\ 
Do 41.9 


-77 13 22 


13.78 


1.47 


0.83 






23 


1 1 


no /ion 

(Jo 4z.y 


-77 10 08 


12.85 


0.95 


0.92 






24 


1 1 
11 


(J8 44.0 


-76 13 29 


12.92 


0.79 


0.56 






25 


1 1 
11 


r\n on A 


-76 34 46 


13.68 


0.95 


1.35 


NIR-14 


2 


26 


11 


r\r\ oo O 

09 62.6 


-76 30 17 


13.00 


0.94 


1.16 






27 


1 1 
11 


09 33.9 


-77 53 33 


13.90 


0.82 


1.19 






28 


1 1 

11 


09 50.7 


-78 04 19 


13.80 


0.77 


1.07 






29 


1 1 
11 


nn o n 

09 52.0 


-76 39 12 


11.67 


0.69 


1.02 


ISO-Chal 217 


1 


30 


1 1 
11 


nn CO A 

09 53.4 


-77 28 36 


12.48 


0.85 


1.24 


ISO-Chal 220 


1 


31 


11 


nn c: y< 1 

09 54.1 


-76 31 11 


12.75 


0.98 


1.11 


ISO-Chal 225 


1 


32 


11 


10 01.6 


-77 47 47 


13.97 


1.23 


1.46 






33 


11 


10 04.5 


-77 48 22 


13.65 


1.21 


1.42 






34 


11 


10 06.2 


-76 40 20 


13.85 


1.06 


1.53 






35 


11 


10 16.8 


-77 43 44 


13.01 


1.38 


1.28 






36 


11 


10 47.6 


-76 32 38 


13.22 


0.84 


0.24 






37 


11 


10 49.4 


-76 38 08 


13.87 


1.19 


0.76 






38 


11 


10 53.2 


-76 16 53 


13.96 


0.78 


0.18 






39 


11 


10 56.9 


-77 13 26 


13.93 


0.77 


1.03 






40 


11 


10 58.2 


-76 17 56 


13.82 


0.77 


0.93 
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Table 3 — Continued 



Star q;(2000.0) 5(2000.0) K H-K J-H Id. Ref. 



41 


11 


11 


08. 


,9 


-76 


49 


12 


13. 


,21 


0. 


,95 


1. 


,54 


42 


11 


11 


21. 


,1 


-78 


05 


20 


13. 


,25 


0. 


,69 


0. 


,64 


43 


11 


11 


27. 


,8 


-76 


50 


17 


13. 


,48 


0. 


,90 


1. 


,09 


44 


11 


11 


35. 


,7 


-76 


53 


26 


13. 


,19 


1. 


,16 


1. 


,79 


45 


11 


11 


42. 


,6 


-78 


05 


21 


13. 


,88 


0. 


,66 


0. 


,53 


46 


11 


11 


48. 


,7 


-76 


12 


21 


13. 


,83 


0. 


,75 


0. 


,49 


47 






^4 


c; 


— 7fi 






X O. 


71 

. ( X 


1 

X . 




n 


Q8 


48 


11 


11 


57. 


,5 


-77 


21 


57 


13. 


,72 


1. 


,38 


1. 


,13 


49 


11 


12 


03. 


,7 


-76 


51 


33 


13. 


,51 


1. 


,27 


1. 


,82 


50 


11 


13 


25. 


,3 


-77 


00 


23 


11. 


,55 


0. 


,98 


0. 


,67 


51 


11 


13 


26. 


,0 


-77 


00 


32 


12. 


,40 


1. 


,20 


0. 


,28 


52 


11 


13 


42. 


,7 


-77 


55 


30 


13. 


,71 


0. 


,45 


0. 


,25 


53 


11 


14 


31. 


,4 


-78 


06 


44 


13. 


,03 


0. 


,98 


0. 


,75 


54 


11 


14 


23. 


,6 


-77 


56 


12 


13. 


,73 


0. 


,84 


0. 


,80 


55 


11 


16 


13. 


,2 


-77 


25 


15 


11. 


,00 


0. 


,85 


0. 


,21 


56 


11 


18 


09. 


,0 


-78 


14 


11 


13. 


,71 


1. 


,03 


0. 


,72 



Note. — Units of right ascension are hours, minutes, and seconds, and 
units of dechnation are degrees, arcminutes, and arcseconds. 

References. — (1) Persi et al. (2000), (2) Oasa et al. (1999) 
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Tablc 4. Candidate Pro-Main Scqiicncc Stars with K > 14.0. 



Star q;(2000.0) 5(2000.0) K H-K J-H Id. Ref. 



57 


10 


40 


47.6 


-77 


36 


49 


15.12 


0.95 


0.73 


58 


11 


04 


17.1 


-77 


25 


49 


15.34 


0.86 


0.65 


59 


11 


04 


27.0 


-77 


28 


02 


14.87 


1.14 


1.02 


60 


11 


04 


30.5 


-77 


34 


41 


14.33 


0.96 


0.89 


61 


11 


04 


32.0 


-77 


54 


09 


15.27 


1.02 


0.88 


62 


11 


04 


35.8 


-77 


29 


09 


15.34 


0.99 


0.94 


63 


11 


04 


36.2 


-77 


50 


12 


14.76 


1.21 


1.26 


64 


11 


04 


38.7 


-77 


24 


21 


14.49 


1.27 


0.87 


65 


11 


04 


39.1 


-77 


40 


05 


14.39 


1.56 


1.18 


66 


11 


04 


48.4 


-77 


49 


20 


15.45 


1.69 


0.85 


67 


11 


05 


37.0 


-76 


51 


18 


15.33 


1.48 


1.04 


68 


11 


05 


41.2 


-76 


38 


38 


15.12 


0.91 


0.58 


69 


11 


05 


42.6 


-77 


26 


30 


14.33 


2.09 


2.67 


70 


11 


05 


43.0 


-77 


31 


58 


14.43 


1.46 


1.16 


71 


11 


06 


33.5 


-77 


52 


26 


14.79 


0.98 


0.74 


72 


11 


06 


35.3 


-77 


21 


10 


15.94 


1.20 


0.95 


73 


11 


06 


49.4 


-77 


34 


37 


14.99 


1.51 


1.78 


74 


11 


06 


51.0 


-77 


11 


33 


15.08 


1.16 


0.64 


75 


11 


06 


54.5 


-77 


25 


54 


15.30 


1.36 


0.70 


76 


11 


06 


59.6 


-77 


18 


29 


14.90 


1.59 


1.52 
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Table 4 — Continued 



Star 0(2000.0) ^KSOOO.O) K H-K J - H Id. Ref. 



11 


ii 


(Y? 
Ui 


c;o 

oy.u 


7Q 
— / 


00 
uu 


io 


1 K OCC 

iO.UO 


1/11 

i.4i 


Q/l 

U.o4 


io 


ii 


Uo 


Zl .6 


— it 


A 7 


oo 
zU 


1 A 7c; 
i4. 1 


1 A A 
i.44 


OQ 

u.yo 


/ y 


ii 


Uo 


Oi.D 


77 
— II 


00 
Uz 


QO 

oU 


1 A QO 
i4.0Z 


1 7C^ 

i. /O 


u.do 


oU 


ii 


Uo 


oi.o 


77 


OR 
UD 


A A 
44 


1 1 
iO.iU 


Q7 
U.o/ 


fi/l 
U.D4 




ii 


no 

uy 


OQ Q 

Uo.o 


7Q 
— / 


1 A 
i4 


zo 


1 c: 00 

io.uy 


1 Q1 

i.oi 


KO 

u.by 


CO 

oz 


ii 


no 

uy 


io.o 


77 
— II 


QO 

cjy 


Q7 


1 c; 1 
iO.iU 


i.ZD 


110 

i.iU 


00 


ii 


no 

uy 


O/l 1 
Z4.i 


7R 
— / D 


Q/l 

o4 


00 


1 /I Q7 

i4.o / 


1 /I 
i.4D 


1 

i.OD 


Q/l 

o4 


ii 


no 

uy 


O/l 7 

z4. / 


7Q 
— / 


io 


c;/i 
04 


1 c: CO 

io.oy 


Qf; 
U.oD 


7Q 

U. / 


Q K 
00 


ii 


no 

uy 


oz. / 


77 
— II 


/I 

4y 


Ofi 

Uo 


1 /I S7 

i4.o / 


1 fiO 

i.oU 


1 Q7 
i.O 1 


OD 


ii 


1 

iU 


01 7 
Ui. / 


77 
— II 


1 
iy 


00 
zU 


1 c: Q c 
iO.oO 


1 oc^ 
i.UO 


1 00 

i.UU 


87 


11 


10 


04.3 


-77 


27 


15 


15.34 


2.01 


0.66 


88 


11 


10 


17.7 


-77 


22 


18 


15.69 


1.31 


1.03 


89 


11 


10 


19.6 


-77 


11 


20 


14.57 


0.94 


0.89 


90 


11 


10 


24.6 


-77 


23 


06 


15.58 


1.93 


0.73 


91 


11 


10 


26.2 


-76 


44 


09 


15.99 


1.18 


0.60 


92 


11 


10 


29.4 


-76 


31 


12 


15.36 


1.37 


0.58 


93 


11 


10 


29.8 


-77 


21 


12 


15.11 


1.43 


0.64 


94 


11 


10 


37.9 


-77 


01 


39 


15.29 


1.21 


0.75 


95 


11 


10 


39.2 


-77 


02 


48 


15.44 


1.45 


1.10 


96 


11 


10 


46.1 


-77 


17 


07 


14.47 


0.99 


0.88 



NIR-11 2 
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Table 4 — Continued 



Star n(2000.0) ^^2000.0) K H-K J - H Id. Ref. 



y / 


ii 


ii 


UD.o 


— io 


nQ 
Do 


K1 


1/17/1 
i4. / 4 


1 nn 
i.UU 


n OQ 
U.yo 


OQ 


i i 


i i 


lie; 


IP. 


A A 
44 


O / 


1 no 
io.Uz 


1 0/1 
i.Z4 


1 OI 

i.zi 


no 

yy 


ii 


ii 


1 O Q 
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Table 4 — Continued 



Star a(2000.0) 5(2000.0) K H-K J-H Id. Ref. 



117 11 17 36.8 -77 45 35 15.18 1.03 0.68 

118 11 18 02.2 -77 41 30 15.61 0.86 0.58 



Note. — Units of right ascension are hours, minutes, and seconds, and 

units of decUnation are degrees, arcminutes, and arcseconds. 
^Star with fuzzy image at K. 



References. — (1) Persi et al. (2000), (2) Oasa et al. (1999) 
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